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Abstract 

The ground-based Atmospheric Radiation Measurement Program (ARM) and NASA 
Aerosol Robotic Network (AERONET) routinely monitor clouds using zenith radiances 
at visible and near-infrared wavelengths. Using the transmittance calculated from such 
5 measurements, we have developed a new retrieval method for cloud effective droplet 
size and conducted extensive tests for non-precipitating liquid water clouds. The under- 
lying principle is to combine a water-absorbing wavelength (i.e. 1640nm) with a non- 
water-absorbing wavelength for acquiring information on cloud droplet size and optical 
depth. For simulated stratocumulus clouds with liquid water path less than 300 gm" 2 
10 and horizontal resolution of 201 m, the retrieval method underestimates the mean effec- 
tive radius by 0.8 pm, with a root-mean-squared error of 1 .7 pm and a relative deviation 
of 13%. For actual observations with a liquid water path less than 450 gm -2 at the 
ARM Oklahoma site during 2007-2008, our 1 .5 min-averaged retrievals are generally 
larger by around 1 pm than those from combined ground-based cloud radar and mi- 
15 crowave radiometer at a 5 min temporal resolution. We also compared our retrievals to 
those from combined shortwave flux and microwave observations for relatively homo- 
geneous clouds, showing that the bias between these two retrieval sets is negligible, 
but the error of 2.6 pm and the relative deviation of 22 % are larger than those found in 
our simulation case. Finally, the transmittance-based cloud effective droplet radii agree 
20 to better than 1 1 % with satellite observations and have a negative bias of 1 pm. Over- 
all, the retrieval method provides reasonable cloud effective radius estimates, which 
can enhance the cloud products of both ARM and AERONET. 


1 Introduction 

Cloud droplet effective radius is one of the most fundamental cloud properties for un- 
25 derstanding cloud formation, dissipation and interactions with aerosol and drizzle (Al- 
brecht, 1 989; Wood, 2000; McComiskey et al., 2009; Kubar et al., 2009). Cloud droplet 
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size is also a crucial determinant of cloud feedback processes, and of the Earth’s ra- 
diative and water energy balance (Slingo, 1990; Wielicki et al. 1995; Stephens, 1999; 
Stephens et al., 2005). While tremendous efforts have been made in providing routine 
cloud droplet effective radii from satellite passive and active measurements (Nakajima 
5 and King, 1990; Han et al., 1994; Kawamoto et al., 2001 ; Chang and Li, 2002; Platnick 
et al., 2003; Roebeling et al., 2006; Mace et al., 2009; Minnis et al., 2011; and many 
others), ground-based retrievals are limited and the discrepancies among different re- 
trievals remain unresolved (Feingold et al., 2006; Schofield et al., 2007; Zhao et al., 
2012 ). 

io Ground-based retrievals of cloud droplet size for liquid clouds are available from 
the Atmospheric Radiation Measurement (ARM) program in Oklahoma, Alaska, and 
the tropical Western Pacific sites (Stokes and Schwartz, 1994). Some retrieval meth- 
ods heavily rely on cloud radar and microwave radiometer measurements (Liao and 
Sassen, 1994; Frisch et al., 1995, 1998; Dong et al., 1998; Dong and Mace, 2003; 
is Wang and Sassen, 2002; Wang et al., 2004), while others utilize a synergy of passive 
radiation measurements. Turner (2007) combined infrared and microwave radiome- 
ter measurements to retrieve cloud droplet radii; their method is limited to cases with 
liquid water paths less than 60 gm " 2 due to saturation of the infrared observations. 
Min et al. (2003) used a least-squared error minimization technique to simultaneously 
20 retrieve cloud optical depth and effective radius using shortwave flux and microwave 
observations. Similarly, Feingold et al. (2006) and Kim et al. (2003) estimated cloud 
droplet radius using liquid water path derived from microwave radiometer and cloud 
optical depth from shortwave flux measurements. 

A number of studies proposed to retrieve cloud droplet size using zenith radi- 
25 ance measurements at visible and near-infrared wavelengths from single instruments 
(Kikuchi et al., 2006; Schofield et al., 2007; Pandithurai et al., 2009; McBride et al., 
2011). The underlying principle for these transmittance-based methods is to com- 
bine a liquid-water-absorbing wavelength with a non-water-absorbing wavelength to 
acquire information on cloud droplet size and optical depth. Kikuchi et al. (2006) 


g 
< /) 
O 
c 

C/5 

C/5 

O’ 


TJ 


0) 

~o 

CD 


D 


C/5 

O 

c 

C/5 

C/5 

o’ 


“0 


Q) 

"O 

CD 


g 

C/5' 

O 

c 

C/5 

C/5 

O' 


Tl 

0 ) 

"O 

CD 


D 

C/5' 

o 

c 

C/5 

C/5 

o’ 

o 

TJ 

CD 

"O 

CD 


ACPD 

12, 19163-19208, 2012 


Cloud droplet size 
and liquid water path 
retrievals 

J. C. Chiu et al. 


Title Page 


Abstract 


Introduction 


Conclusions 


References 


Tables 


Figures 




Back 


Close 


Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 



19166 





and Pandithurai et al. (2009) used zenith radiance measurements at 1020, 1600 and 
2200 nm wavelengths to simultaneously retrieve cloud optical depth and effective radius 
with a 1 min temporal resolution. For water clouds with effective radii ranging between 
1-6 pm, the retrieved cloud droplet radii from zenith radiance data were generally larger 
5 than those retrieved from cloud radar in Kikuchi et al. (2006), but vice versa in Pandithu- 
rai et al. (2009). Schofield et al. (2007) used zenith radiances at 0.9-1 .7 pm to observe 
a homogeneous cloud over Barrow, Alaska, and the cloud droplet radii derived were 
about 3 pm larger than cloud radar retrievals. Recently, McBride et al. (2011) developed 
a novel method for cloud droplet size retrievals using hyperspectral measurements 
io from the ARM shortwave spectrometer at 1 s resolution; they enhanced the sensitivity 
of transmittance observations to effective radius by using the spectral slope derived 
from measurements between 1 565 and 1 634 nm. The new McBride et al. (201 1 ) spec- 
tral method led to good agreement in both effective radius and liquid water path, when 
compared to satellite retrievals and ground-based microwave observations, 
is Figure 1 summarizes the aforementioned shortwave-transmittance-based retrievals. 
Overall, the comparison of transmittance-based retrievals to cloud radar retrievals is 
less conclusive and depends strongly on radar retrieval methods. On the other hand, 
cloud droplet radii retrieved from ground-based passive observations tend to be smaller 
than those from satellite reflectance. This is because the ground- and satellite-based 
20 retrievals obtain their information from different levels within the cloud layer. Platnick 
(2000) used the number of photon scattering events to weight the contribution from 
each level to the overall size determination, and found that the radiative contribution 
mainly comes from the upper portion of the cloud layer for satellite-based retrievals. 
In contrast, for ground-based retrievals, the contribution to the reported size comes 
25 from all levels within the cloud layer. Models from quasi-adiabatic parcel to Large Eddy, 
as well as in situ measurements, generally show that liquid cloud droplets grow from 
cloud base to cloud top, which would explain why the droplet radii from ground-based 
retrievals are smaller than from satellite-based retrievals. 
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The Aerosol Robotic Network (AERONET) recently introduced a new observation 
strategy, called “cloud mode”, which dramatically increases the number and variety 
of cloud observations on a global scale (Chiu et al., 2010). AERONET is comprised 
of sun/sky radiometers with a 1.2° field-of-view designed for aerosol retrieval. When 
5 clouds completely block the sun - making measurements of aerosol properties less 
practical - the radiometer switches to “cloud mode” and performs 10 zenith radiance 
measurements at 9 s intervals, then goes to sleep for 15 min. Among the sites employ- 
ing cloud mode shown in Fig. 2, Table 1 lists 31 AERONET sites where zenith radiance 
measurements also include the 1640nm wavelength. A combination of the 1640nm 
io liquid water-absorbing wavelength with the others (e.g. 440, 675, 870 or 1 020 nm) can 
be used to retrieve cloud effective radius to potentially enhance the cloud mode product 
that currently provides retrievals only of cloud optical depth. 

This paper presents a new retrieval method for cloud effective droplet radius, us- 
ing zenith radiance measurements from ARM and from the AERONET cloud mode 
15 operation. In Sect. 2, we describe our retrieval method and evaluate its performance 
on a stratocumulus cloud from a large eddy simulation using Monte Carlo radiative 
transfer calculations. In Sect. 3, we introduce several ancillary observational datasets 
and discuss the results of the intercomparison. In particular, intercomparisons involving 
satellite retrievals are placed within the context of Fig. 1. Finally, Sect. 4 summarizes 
20 this paper. 


2 Methodology for retrieving cloud droplet size 
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2.1 The retrieval method 
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For clouds over a Lambertian surface, the ground-based zenith radiance / at a wave- 
length A is a function of cloud optical depth and effective radius. In general, it can be 
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(a — f ^eff> t^O’Px) 


where r c is cloud optical depth; A c is effective cloud fraction; r eff is cloud effective 
radius; /u 0 is the cosine of solar zenith angle; and p x is the albedo of the underlying 
surface. To retrieve three parameters of interest (i.e. r c , A c , r eff ) in Eq. (1), we extended 
the method in Chiu et al. (2010) that used 440 and 870 nm wavelengths by adding 
a 1640nm water-absorbing wavelength. 

Relationships between zenith radiance and cloud optical depth at 870 and 1640nm 
for droplet sizes of 8 and 16 pm are illustrated in Fig. 3. It is known that the larger the 
cloud droplets, the stronger the absorption and forward scattering. This initiates a com- 
peting process; stronger absorption for larger droplet sizes reduces zenith radiance 
reaching the surface, while stronger forward scattering enhances it. As a result, the 
zenith radiance at 870 nm wavelength increases with droplet size due to the dominant 
factor of forward scattering. In contrast, the radiance at 1640 nm wavelength decreases 
with droplet size due to the dominant factor of absorption; however, the sensitivity of the 
zenith radiance to droplet size is reduced by the increase in forward scattering (Rawlins 
and Foot, 1990; Platnick, 2000). This implies that uncertainties in zenith radiance mea- 
surements and surface albedo estimates may have a non-negligible impact on droplet 
size retrievals, and, in turn, need to be accounted for in the retrieval method. 

In our retrieval method, we assumed 5% uncertainty in zenith radiance measure- 
ments at all wavelengths, 5 % uncertainty in surface albedo at 870 and 1 640 nm wave- 
lengths, and 10% uncertainty in surface albedo at 440 nm wavelength. These un- 
certainties, normally distributed and estimated from Holben et al. (1998) and Schaaf 
et al. (2002), are used to perturb the observed zenith radiance and surface albedo 
estimate, resulting in an overall relative input uncertainty of ~17%. We compare the 
perturbed zenith radiance to calculated lookup tables and search for possible solutions. 
The lookup tables were computed from the discrete-ordinate-method radiative transfer 
model (Stamnes et al., 1988) over reasonable ranges of cloud optical depth, effective 

19169 
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cloud fraction, and cloud effective radius ranging from 4 to 20 pm, typical for ARM sites 
over Oklahoma, Alaska and Darwin (Zhao et al., 2012). We also postulate that the so- 
lution is found when zenith radiances agree with the lookup table to within 10% at both 
440 and 870 nm wavelengths, and sort those solutions based on errors in the zenith 
radiance at 1640nm; the smaller the error, the better the solution. Finally, we take the 
mean of the best five solutions to produce the best solution for this particular set of the 
perturbed zenith radiance and surface albedo. 

For each cloud-mode measurement, the above procedure is repeated 40 times using 
randomly generated perturbations. The final retrievals for cloud optical depth and effec- 
tive radius are given as the mean of these 40 repetitions with an uncertainty estimated 
by the standard error. Although the thresholds of five best solutions and 40 repetitions 
are somewhat arbitrary, the mean bias and the root-mean-squared difference change 
insignificantly (about 2%) when the thresholds increase to 10 best solutions and 1000 
repetitions in the following simulation test. 

2.2 Simulations setup and evaluation metrics 

To evaluate the performance of our retrieval method, we tested a stratocumulus cloud 
from the Intercomparison of 3-D Radiative Codes (I3RC; Cahalan et al., 2005). This 
cloud case was based on large eddy simulations of FIRE-1 stratocumuli (Moeng et al., 
1996; Hinkelman et al., 2007) and consisted 64x64x 16 cells. In this simulation test, we 
used two horizontal resolutions to investigate at which scale the retrieval method per- 
forms reasonably well. One set retains the original spatial resolution of the I3RC input 
data, 67x67x25 m. By expanding grid cells horizontally but keeping the same values as 
the original smaller cells, we obtained the other set with resolution of 201 x 201 x 25 m, 
close to the minimum horizontal scale where radiation measurements retain the statis- 
tics of the cloud structures (Marshak et al., 1995). Using the latter resolution, Fig. 4 
shows an image of the overall cloud optical depth, and vertical features for a particular 
transect at 3.1 km in the Y direction. As Fig. 4b shows, cloud geometric thicknesses 
are about 400 m and effective droplet radii range between 4 and 20 pm. 

19170 
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Zenith radiances for this particular transect were simulated using a backward Monte 
Carlo model (similar to the model UMBC2 in Cahalan et al., 2005). We assumed that 
the solar azimuth and zenith angles are respectively 0° and 45°. The surface albedo 
values at 440, 870 and 1640nm were respectively 0.05, 0.3 and 0.25, common for 
5 a vegetated surface based on the Collection 5 products of the Terra and Aqua Moderate 
Resolution Imaging Spectroradiometer (MODIS) combined data set at 500 m resolution 
(Schaaf et al., 2002). 

Using this transect, we retrieve optical depth, cloud effective radius, and liquid water 
path (LWP) using the simulated zenith radiances. The former two parameters can be 
io obtained directly from our method, but LWP needs to be calculated. One approach is 

LWP = |p w Velf> (2) 

where p w is the density of water; this equation is based on the assumption that liquid 
water content (LWC) is constant in the vertical (Stephens, 1978). Another approach is 

is LWP = ^p w T c r eff , (3) 

which assumes that the droplet number concentration is approximately constant and 
the liquid water content increases linearly with height (Wood and Hartmann, 2006). 
Equation (2) is exactly 20 % larger than Eq. (3). Since in the I3RC input data the optical 
20 depth for each level was derived from Eq. (2), we used the same equation to calculate 
LWP from retrieved cloud optical depth and effective radius in this test. 

Now we need to determine the “true” column effective radius. Recall that the I3RC 
input data includes LWP, r c , and r en for each level within the cloud layer. While the true 

total LWP and optical depth can easily be calculated by integrating the entire column, 

25 the “column cloud effective radius” cannot be unambiguously defined. One option is 
to use Eq. (2) to calculate r eff from the column-integrated LWP and the optical depth, 
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which we refer to as r eff C onstLwc- Another option is to use the mean effective radius r eff : 


N 


^eff = 


N 


r, 


eff, level’ 


levels 1 


( 4 ) 


where r eff j eve! is the effective radius at each cloud level given by the I3RC input data and 
5 N is the number of the cloud levels in the cloud layer. This estimation represents the 
fact that for ground-based retrievals the contribution to the overall size determination is 
relatively uniform over the cloud layer (see Sect. 1). 

Once the truth and the retrievals are generated, we use several evaluation metrics 
to characterize the agreement between the two: the mean bias (BIAS; the difference 
io between the means), the root-mean-squared difference (RMSD), the linear correlation 
coefficient (CORR) and the relative deviation (RD, in %) expressed as the ratio of the 
average absolute deviation to the mean value, i.e. 


M 

2 | ^eff, retrieved,/' — ^eff.true,/ 1 

RD = — x100, (5) 

M 

^ ^eff.tru e,i 
i = 1 


where M is the total number of data points of interest; /"eft, true,/ ancl ^.retrieved,/ are the 
15 true and retrieved effective radius at a given data point /, respectively. 
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2.3 Evaluation results from simulations 
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Let us first focus on the overall performance of the retrievals from simulations at 201 m 
horizontal resolution. Table 2 and Fig. 5 show that the liquid water path retrievals and 

the truth have similar means, standard deviations and distribution percentiles; the cor- 

—2 —2 

20 responding mean bias and RMSD are -5 g m and 1 9 g m , respectively. The RMSD 
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is comparable to the uncertainty of 20-30 gm" 2 in liquid water path retrievals from two- 
channel microwave radiometer measurements (Marchand et al., 2003; Crewell and 
Lohnert, 2003). Similarly, cloud optical depth retrievals agree well with the truth; the 
corresponding mean bias and RMSD are 0.7 and 2, respectively. 

5 Retrieval errors in cloud effective radius depend on which definition of the true col- 
umn effective radius is used, i.e. r e ff, constLwc from Eq. (2) or 7 eff from Eq. (4). Table 2 
and Fig. 6 show that the retrieval errors are smaller and the box plots agree better 
when we use r eff as the truth instead of /" e ff, constLwc- However, for both situations, the 
correlations between the retrieved and true effective radii are weak for reasons given 
io in Appendix A. 

To check retrievals in more detail, Fig. 7 shows the box plots for all locations along 
the transect at 3.1 km shown in Fig. 4. Each box represents the statistical distribution 
of the effective radii for all levels at a given location in the x-direction. For 62 cloudy 
points and two clear-sky points, only one cloudy point is not retrievable. Among the 
15 retrievable cloudy points, 80% fall between the 25th and 75th percentiles, and 60% 
are lower than the median. Overall, the majority of retrievals agree well with the true 
effective radius distributions statistically, although it remains unclear which definition of 
column effective radius is better. 

Finally, to evaluate how retrieval performance changes with horizontal resolution, 
20 error statistics for retrievals at the 67 m horizontal resolution are summarized in Ta- 
ble 2. Compared to results at the 201 m horizontal resolution, the relative deviation in 
retrieved LWP at the smaller scale increases by 10%, while the relative deviations in 
retrieved cloud optical depth and effective radius increase by 5 %. 

To summarize, a better agreement has been achieved at the 201 m horizontal scale, 
25 where retrieved cloud effective radius has uncertainty 15% and RMSD 1.4-1 .7 pm, 
while retrieved cloud optical depth has uncertainty 15% and RMSD 2. 
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3 Comparison of cloud-mode retrievals to ARM and MODIS products 

For real world applications, we use 1 .5 min-averaged zenith radiances from cloud-mode 
measurements to retrieve cloud optical depth and effective radius. Surface albedo es- 
timates are obtained from MODIS collection 5 products (Schaaf et al., 2002). When 
5 these are not available, we use an optimized climatological database derived from 
2000-2004, 16 day average MODIS surface reflectivity product at 1 min resolution 
(Moody et al., 2005, 2007; Eck et al., 2008). The resulting cloud-mode retrievals are 
then compared to several products from ARM and MODIS during May 2007-June 
2008. We focus on the ARM Oklahoma site due to the superior availability of its data 
io products, including LWP from microwave radiometer measurements, and cloud effec- 
tive radius from flux and radar measurements. 

3.1 LWP comparison to ARM microwave- and infrared-based retrievals 

In this section, we evaluate cloud-mode LWP retrievals against the ARM Archive MWR- 
RET product (Turner et al., 2007). MWRRET optimizes LWP estimates by blending 
is radiative transfer calculations with radiosonde measurements and with microwave ra- 
diometer measurements having a 5.9° field-of-view and 20 s time resolution. This 
physical-based method significantly reduced the clear-sky bias, a common problem 
in all earlier retrieval techniques. In general, the retrieval uncertainty from MWRRET is 
20-30 gm“ 2 . 

20 Since microwave radiometers mainly detect liquid clouds and retrievals are unreli- 
able during precipitating periods, we excluded rainy cases using rain gauge measure- 
ments, and ice cloud cases using Cloudnet products (Illingworth et al., 2007). Cloud- 
net provides ice water content at 30 s time resolution and 90 m vertical resolution, us- 
ing cloud radar reflectivity measurements and temperature profile information (Hogan 
25 et al., 2006). Although the retrieval uncertainty in ice water content is 25-50 % (Heyms- 
field et al., 2008), a threshold of zero ice water content has been applied for our ice 
cloud exclusions. 
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The LWP values from the MWRRET product during May 2007-June 2008 range 
between 10 and 450 gm" 2 with a mean of 117grrT 2 and one standard deviation of 
83 gm“ 2 . For intercomparison, we used Eqs. (2) and (3) to calculate retrieved LWP. 
Table 3 shows that the use of Eq. (3) significantly reduces the bias from 28 down to 4 
5 g rrf 2 , and slightly reduces the RMSD from 60 to 50 gm -2 . This indicates that cloud 
water content more likely increases with height at the ARM Oklahoma site, consistent 
with the findings of McBride et al. (2011). 

To take a more detailed look, Fig. 8a shows that LWP retrievals calculated from 
Eq. (3) correlate well to MWRRET. Figure 8b shows that histograms from the two 
io methods peak at 50-1 00 gm -2 , but the occurrence of low LWP values is much more 
frequent in MWRRET than in cloud mode. Although this discrepancy might be partly 
due to the difference in field-of-view and observation strategy, it requires other datasets 
to understand the cause. 

3.2 Cloud effective radius comparison to ARM radar-based retrievals 

is A continuous dataset for cloud effective radius at a time resolution of 5 min and a verti- 
cal resolution of 90 m is available in the ARM Archive, based on merging measurements 
of cloud radar, microwave radiometer and radiosonde soundings (Mace et al., 2006). 
For daytime retrievals, the effective radius of liquid clouds at cloud base is estimated 
from an empirical relationship using liquid water path and solar transmission (Dong 
20 et al., 1998). The effective radius at cloud base is then used to estimate the vertical 
distribution of effective radius based on the observed radar reflectivity profile (Dong 
and Mace, 2003). We call this dataset “ARM Mace” hereafter. 

Cloud radar data in ARM Mace were excluded from comparison based on several 
factors: (1) when solar zenith angle was greater than 70°; (2) when rain gauge data 
25 showed non-zero rain rate; (3) when the window of the microwave radiometer was wet; 
and (4) when ice water contents were non-zero. After these exclusions, mean effective 
radii were calculated by averaging over all levels within the cloud layer. 
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Occurrence histograms of cloud effective radii from cloud mode and coincident ARM 
Mace retrievals are shown in Fig. 9. The occurrence frequency of ARM Mace 5 min 
average effective radii peaks at 6-8 pm, while the frequency of cloud-mode retrievals 
peaks at 8-1 Opm. Cloud-mode retrievals also show a much higher occurrence fre- 
quency at 10-1 4 pm compared to the Mace dataset. This in turn leads to a significant 
difference in the overall mean in Fig. 9b, where we see that effective radii retrieved from 
the ARM Mace and cloud-mode datasets are respectively (7.6 ±1 .8) and (8.5 ±2.2) pm. 
The bias of ~ 1 pm is smaller than those reported in Schofield et al. (2007) and Pan- 
dithurai et al. (2009), but similar to those reported in Kikuchi et al. (2006), as discussed 
in Sect. 1. 

Since cloud properties are scale-dependent, the bias of ~ 1 pm is partly due to the 
fact that these two retrievals represent different time scales. Unfortunately, 5 min aver- 
age cloud-mode retrievals cannot be derived due to less frequent cloud-mode observa- 
tions, which are dependent on the ~ 1 5 min direct sun measurement frequency (Holben 
et al., 1998). To identify the source of the discrepancy, we will need to either change 
cloud-mode observation strategy or do radar retrievals with higher time resolution. 

3.3 Cloud effective radius comparison to ARM flux-based retrievals 

In this section, we evaluate cloud-mode effective radii against the ARM product that 
is based on shortwave flux measurements. The ARM multifilter rotating shadowband 
radiometer (MFRSR), with a hemispheric FOV, provides 20 s averages of both direct 
and diffuse solar flux in narrow bands centered at 41 5, 500, 61 5, 673, 870, and 940 nm. 
Min and Harrison (1996a) used direct and diffuse transmittance at 415nm to estimate 
cloud optical depth with an initial default effective radius of 8 pm. With additional liquid 
water paths retrieved from the ARM microwave radiometers, cloud effective radius and 
optical depth can be simultaneously retrieved by minimizing least-squares errors in 
radiance along with an adjoint radiative transfer method (Min and Harrison, 1996b; Min 
et al., 2003). If the error minimum can’t be found, the default 8-pm size is reported 
and is excluded in our intercomparison. In addition, to match the temporal resolution 
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of AERONET cloud-mode retrievals, effective radius values from the ARM dataset are 
averaged to a 1.5 min time window. We call this dataset “ARM Min” hereafter. 

Results for a 15 June 2007 case are investigated as shown in Fig. 10. The plot 
of lidar backscatter coefficients shows that liquid layers severely attenuate the signal 
5 during 1 4:00-1 7:00 and 1 8:00-1 9:00 UTC, indicating the presence of optically thick low 
clouds. It also shows occasional high clouds during 17:00-18:00 and after 19:00 UTC. 
Figure 10b and c show that in the time periods of continuous low clouds, retrievals of 
both cloud effective radius and optical depth from cloud mode and ARM Min datasets 
agree reasonably well, except near 17:00 UTC when cloud optical depth decreased 
io dramatically from 100 to 25 and cloud-mode effective radius is higher by ~4gm. For 
time periods when clouds are less homogeneous (like the case around 1 7:00 UTC just 
mentioned), droplet radii from the two datasets are significantly different. 

To extend our comparison to all non-precipitating and relatively homogeneous liquid 
clouds, we first used rain gauge data to exclude rainy periods. We further excluded 
15 the time period when cloud gaps were present, by requiring the ARM Min retrievals of 
optical depth to be continuously greater than 1 0 for at least one hour. We also excluded 
time periods when cloud optical depths fluctuated significantly, by checking whether 
the corresponding standard deviation exceeded 10. Finally, we excluded overlying ice 
cloud cases using Cloudnet products and the procedure described in Sect. 3.1. 

20 For non-precipitating and relatively homogeneous liquid clouds, Fig. 11a shows that 
cloud optical depths agree well with those from the ARM Min dataset, with a correla- 
tion coefficient 0.96. Generally, cloud optical depths from cloud mode are larger than 
those from ARM Min by a mean difference of 3; the corresponding RMSD and rela- 
tive deviation are 5.5 and 13%, respectively (Table 4). In contrast, the correlation in 
25 cloud effective radius shown in Fig. 11b is low (0.3), even though the mean and stan- 
dard deviation of effective radius from the two are nearly identical. This is explained in 
Appendix A. The corresponding RMSD of 2.6 pm and the relative deviation of 22% in 
effective radius retrievals are larger than the uncertainty estimated from our simulation 
test (Table 2). 
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3.4 Cloud effective radius comparison to MODIS reflectance-based retrievals 

In this section, we evaluate cloud-mode effective radii against MODIS Level 2 cloud 
products (Collection 5), which provide cloud effective droplet radius and cloud phase 
at 1 km resolution. We used droplet radii retrieved from the MODIS 2.1 pm wavelength, 
5 which agreed with those retrieved from the 3.7 pm wavelength to within 2 pm for rela- 
tively homogeneous clouds (Zhang and Platnick, 201 1 ). For our intercomparison, rela- 
tively homogenous cases were selected when cloud optical depths from the ARM Min 
dataset were continuously greater than 10 for at least one hour. 

Similar to Chiu et al. (2010), we have chosen 1 h time windows centered at the 
io MODIS overpass times for our intercomparison; the number of MODIS pixels was de- 
termined by the wind speed at the cloud base height. In addition, we used MODIS re- 
trievals only when their corresponding cloud phase is liquid water. We also used cloud- 
mode retrievals only when MODIS measurement times were ice-free based on Cloud- 
net products. As a result, eight Terra and Aqua overpasses during May-December 
is 2007 are used and listed in Table 5. 

The box plot in Fig. 12 shows that all cloud-mode points fall into the whiskers if not 
always the boxes of MODIS retrievals, and 75% of points are in the lower percentile 
of the satellite retrieval distributions. As discussed in Sect. 1, satellite-based effec- 
tive radius retrievals are mainly determined by droplet sizes nearer to the cloud tops, 
20 while ground-based retrievals are affected by all cloud layers. Thus, we expect cloud- 
mode retrievals to be distributed in the lower 50th percentile of the satellite retrieval 
range, which is consistent with the behavior seen in Fig. 12. Overall, Table 6 shows 
a bias of -0.9 pm and a relative deviation of 1 1 % between the two sets of retrievals. 
Note that uncertainties in cloud-mode retrievals can be 15% based on our simulation 
25 test, and uncertainties in MODIS retrievals can be 5-20% due to cloud inhomogene- 
ity (Nakajima et al., 1991; Platnick and Valero, 1995). Painemal and Zuidema (2011) 
also reported a 15-20% overestimate of MODIS retrievals compared to in-situ data 
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for stratocumulus clouds. These reported uncertainties help conclude that cloud-mode 
retrievals agree reasonably well with MODIS for relatively homogenous clouds. 

Finally, the means of cloud-mode and MODIS retrievals for each overpass are plotted 
in Fig. 1 3, along with those reported in literature (similar to Fig. 1 ). In general, our over- 
cast cases fall into the radius range between 9 and 13 pm. Some cloud-mode points 
are close to the one-to-one diagonal line, similar to the points in McBride et al. (2011); 
some have slight deviations from MODIS retrievals, similar to the point in Feingold 
et al. (2006). 


4 Summary 
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We developed a new retrieval method for cloud effective droplet radius using ground- 
based zenith radiances taken from the ARM program and the NASA AERONET. The 
method takes advantage of the additional 1640nm wavelength that has recently be- 
come standard for monitoring aerosol properties in AERONET. Specifically, we used 
440, 870 and 1640nm wavelengths to simultaneously retrieve cloud optical depth and 
effective droplet size, based on the fact that radiance at the first two wavelengths is 
mainly determined by cloud optical depth, while at the last wavelength it is mainly de- 
termined by droplet size. 

We conducted extensive tests using simulations, ground- and satellite-based mea- 
surements for non-precipitating liquid water clouds. First, we tested the retrieval method 
on a simulated stratocumulus cloud at two different spatial scales, 67 and 201 m. The 
transect we took for evaluation represented clouds with LWP less than 300 gm" 2 , and 
effective radii ranging between 4 and 20 pm. Overall, a better agreement was achieved 
for the 201 m spatial scale; the retrieval uncertainty in cloud effective radius is 13% and 
the RMSD is 1.7 pm. At the 67 m spatial scale, errors of LWP increase by 10%, and 
errors of cloud effective radius and optical depth increase by 5 %. 

Second, we compared our retrievals to those from cloud radar and microwave mea- 
surements. For clouds with liquid water path less than 450 gm -2 , better agreement has 
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been achieved under the assumption that liquid water content increases linearly with 
height; the retrieval uncertainty in liquid water path is 30% and the root-mean-squared 
error is 50 gm" 2 . In addition, coincident cloud-mode and ARM cloud radar retrievals 
have shown that the latter has a narrower peak at 6-8 pm than the former which 
5 peaks at 6-1 0 pm. The overall bias is 1 pm, comparable to that reported in Kikuchi 
et al. (2006). 

Third, we compared our retrievals to those from flux and microwave combined mea- 
surements for relatively homogeneous clouds. Our retrieval uncertainty in liquid water 
path is 26 % and the RMSD is 46 g m -2 , similar to the finding for all-liquid-water clouds. 
io This is not surprising because effective radii from the flux-based method were con- 
strained by microwave observations. The retrieval uncertainty in cloud effective radius 
is 22 % and the RMSD is 2.6 pm; these values are larger than those suggested by the 
simulation case. Interestingly, even though a low correlation coefficient is obtained in 
cloud effective radius between the two retrieval sets, the bias of 0.1 pm between the 
is two retrieval sets is negligible. This is explained in Appendix A. 

Finally, we compared our retrievals to satellite-based retrievals that have been widely 
used for cloud study. Cloud-mode retrieval uncertainty in effective droplet radius is 1 1 % 
and the RMSD is 1.8 pm. In particular, cloud-mode retrievals have a negative bias of 
0.9 pm, compared to satellite data. This negative bias can be attributed to the fact 
20 that cloud droplet effective radius increases with height for liquid layers under quasi- 
adiabatic conditions; ground-based retrievals receive an equal contribution to the layer- 
averaged cloud droplet effective radius from all levels within the cloud layer, whereas 
the dominant contribution for satellite-based retrievals is from the cloud top. 

Overall, for non-precipitating liquid water clouds, the retrieval uncertainty is 11-22% 
25 for cloud effective radius and 14-30% for liquid water path based on a simulated stra- 
tocumulus case and on ground and satellite measurements. For the ARM Oklahoma 
site, the assumption of a linear increase in liquid water content with height appears 
to be suitable for all clouds. At other sites, caution must be exercised, with additional 


g 

C/3 

o 

c 

C/5 

C/5 

O’ 


TJ 

03 

■o 

CD 


U 


C/5 

o 
c 
< / > 
C/5 

o’ 


“0 


03 

"O 

CD 


U 

(/)' 

o 

c 

C/5 

C/5 

O' 


Tl 

03 

“O 

CD 


D 

C/)' 

o 

c 

C/5 
c n 
o’ 


"D 

03 

■o 

CD 


ACPD 

12, 19163-19208, 2012 


Cloud droplet size 
and liquid water path 
retrievals 

J. C. Chiu et al. 


Title Page 


Abstract 


Introduction 


Conclusions 


References 


Tables 


Figures 




Back 


Close 


Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 



19180 





intercomparisons against microwave observations necessary to determine an appro- 
priate assumption for the vertical distribution of liquid water content. 


Appendix A 


The correlation between the truth and the retrievals is determined not only by the re- 
5 trieval accuracy, but also by the dynamic range of the retrieval variable. In our simulation 
test, the majority of cloud effective radii range between 8 and 12 pm. For such small 
dynamic range, achievement of a high correlation coefficient requires high retrieval ac- 
curacy, which is illustrated next using synthetic data. 

Suppose that retrievals agree to the truth within an error of e. We randomly generated 
io retrievals by 

/"retrieve = 'true + £-{2a-A), (A1) 

where a is a random variable uniformly distributed between 0 and 1; r true andr^i^g 
represent the truth and its retrieved counterpart, respectively. Figure A1 shows that 
is the correlation coefficient between r true and /"retrieve decreases with increasing error of 
e. When the error is comparable to the dynamic range, the correlation automatically 
becomes weak as shown in Fig. Ale and d. Our retrieval uncertainty is on the order 
of 1.5 pm and comparable to an overall uncertainty in input data of 17%. The error 
interval (~ 3 pm) however is comparable to the dynamic range (~ 4 pm) of the effective 
20 radius, forcing the correlation coefficient to always be low. In other words, retrievals can 
be accurate (i.e. comparable to input uncertainties) but weakly correlate to the truth. 
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Table 1 . A list of AERONET cloud-mode sites where zenith radiance measurements at 1 640 nm 
are available in 201 1 . 


Site Country Longitude (°) Latitude (°) 


130.9 - 12.4 

149.1 - 35.3 


ARM Darwin 

Canberra 

Ragged_Point 

XiangHe 

Camaguey 

Carpentras 

Lille 

Palaiseau 

Paris 

REUNION_ST_DENIS 

Hamburg 

ARM_Gan_lsland 

ARMJMainital 

IMAA_Potenza 

Nauru (ARM) 

Manus (ARM) 

Cabo_da_Roca 

Graciosa (ARM) 

Autilla 

Burjassot 

La_Laguna 

EPA-NCU 

Barrow 

Brookhaven 

Cart_Site (ARM) 

GSFC 

Monterey 

SERC 

UMBC 

Univ_of_Houston 

NGHIA_DO 


Australia 

Australia 

Barbados 

China 

Cuba 

France 

France 

France 

France 

France 

Germany 

India 

India 

Italy 

Nauru 

Papua New Guinea 

Portugal 

Portugal 

Spain 

Spain 

Spain 

Taiwan 

United States 

United States 

United States 

United States 

United States 

United States 

United States 

United States 

Vietnam 


- 59.4 

13.2 

117.0 

39.8 

- 77.8 

21.4 

5.1 

44.1 

3.1 

50.6 

2.2 

48.7 

2.3 

48.9 

55.5 

- 20.9 

10.0 

53.6 

73.1 

- 0.7 

79.5 

29.4 

15.7 

40.6 

166.9 

- 0.5 

147.4 

- 2.1 

- 9.5 

38.8 

- 28.0 

39.1 

- 4.6 

42.0 

- 0.4 

39.5 

- 16.3 

28.5 

121.2 

25.0 

- 156.7 

71.3 

- 72.9 

40.9 

- 97.5 

36.6 

- 76.8 

39.0 

- 121.9 

36.6 

- 76.5 

38.9 

- 76.7 

39.3 

- 95.3 

29.7 

105.8 

21.0 
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Table 2. Statistics of true and retrieved liquid water path (LWP), cloud optical depth t c and 
effective radius (^ e ff,constLwc ar) d ^eff) from the simulations at 201 m and 67 m horizontal resolution 
presented in Sect. 2; mean and standard deviation are included in each parenthesis. The other 
columns represent mean bias (BIAS), root-mean-squared difference (RMSD), relative deviation 
(RD in %) and correlation coefficient (CORR). 


Variable 

Truth 

Retrieval 

BIAS 

RMSD 

RD (%) 

CORR 

Based on simulations at 201 m resolution 

LWP (g irf 2 ) 

(81 ±52) 

(76 ± 42 ) a 

-5 

19 

14 

0.95 

T c 

(10.7 ± 5.8) 

(11. 4 ±6.0) 

0.7 

2.0 

14 

0.95 

^eff.constLWC (M^) 

(10.7 ± 1 .6) b 

(9.9 ±1.2) 

-0.8 

1.7 

13 

0.43 

'eft (t^ m ) 

(9.91 ±1.25) c 

(9.85 ±1.19) 

-0.06 

1.4 

12 

0.32 

Based on simulations at 67 m resolution 

LWP (g irf 2 ) 

(81 ±52) 

(68 ± 30) a 

-13 

34 

25 

0.86 


(10.7 ± 5.8) 

(11. 2 ±5.7) 

0.5 

3.0 

19 

0.88 

Teff.constLWC (M*Tl) 

(10.7 ± 1 .6) b 

(9.6 ±1.6) 

-1.1 

2.7 

20 

0.20 

h tt (Mm) 

(9.9 ± 1 .3) c 

(9.6 ±1.6) 

-0.3 

2.4 

18 

0.30 


a LWP values are calculated using retrieved r c and effective radius assuming constant liquid water content in 
the vertical. 

b Calculated from the column-integrated LWP and t c assuming constant liquid water content in the vertical. 
c Calculated from averaging effective radii over the cloud levels. 
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Table 3. Statistics for comparison of liquid water path (gm -2 ) between the new cloud-mode 
retrievals and retrievals from the ARM microwave radiometer (MWR) during May 2007-June 
2008. Descriptions of these metrics can be found in Table 2 (but here BIAS is computed treating 
the MWR as truth). The ARM product used for this analysis is called MWRRET in the archive. 
Cloud-mode LWP is calculated from retrieved cloud optical depth and effective radius using two 
approaches. 


ARM MWR LWP 

Cloud-mode LWP 

BIAS 

RMSD 

RD 

CORR 

(gm- 2 ) 

(gm' 2 ) 

(gm -2 ) 

(gm- 2 ) 

(%) 


(1 17 ± 83) 

(145 ± 82) a 

28 

60 

38 

0.8 


(121 ±68) b 

4 

50 

30 

0.8 


a Calculated from Eq. (2), assuming constant liquid water content in the vertical. 
b Calculated from Eq. (3), assuming an approximately constant cloud droplet concentration 
and a linear increase in liquid water content with height. 
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Table 4. Statistics for comparison of liquid water path, cloud optical depth and effective radius 
between the new cloud-mode retrievals and the ARM flux-based retrievals. Cases are limited 
to non-precipitating and relatively homogeneous clouds during May 2007-June 2008. 


Variable 

ARM Min 

Cloud mode 

BIAS 

RMSD 

RD (%) 

CORR 

LWP (g m -2 ) 

(132 ±54) 

(148 ±71)* 

16 

46 

26 

0.80 

Cloud optical depth 

(30 ±13) 

(33 ±15) 

3 

5.5 

13 

0.96 

Effective radius (pm) 

(8.2 ±2.2) 

(8.3 ±2.2) 

0.1 

2.6 

22 

0.30 


* Calculated from the retrieved cloud optical depth and effective radius assuming an approximately constant 
cloud droplet concentration and a linear increase in liquid water content with height. 
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Table 5. Information on Terra/Aqua MODIS overpasses in 2007 used for evaluations of cloud- 

n 

mode retrievals. 
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* Estimated from ARM’s Active Remotely Sensed Clouds Locations (ARSCL) product. 
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Table 6. Statistics for comparison of cloud effective droplet radius (pm) to MODIS retrievals 
using eight Terra and Aqua overpasses in May-December 2007. The metrics include the mean 
bias (BIAS), root-mean-squared difference (RMSD), relative deviation (RD in %) and correlation 
(CORR). 


MODIS (pm) 

Cloud mode (pm) 

BIAS (pm) 

RMSD (pm) 

RD (%) CORR 

(1 1 .2 ± 1 .5) 

(10.3 ±1.6) 

-0.9 

1.8 

11 0.5 
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Fig. 1- Retrieval intercomparison of cloud effective droplet radii reported in selected literature. 
Data points using cloud radar retrievals as the source effective radii are plotted with dots, while 
data using MODIS retrievals as the source are plotted with triangles. Cloud effective radii from 
Feingold et al. (2006) were retrieved from combined flux and microwave measurements, and 
represented 3 h statistics for stratus clouds at the ARM Oklahoma site in 2003. Retrievals from 
Pandithurai et al. (2009) were based on multichannel zenith radiance and represented 3-day 
statistics for low-level, overcast ice-free water clouds at Okinawa Island, Japan in 2008; the use 
of liquid water-absorbing wavelength (i.e. 1 .6 or 2.2 pm) slightly changed the retrieved cloud 
effective radii. Similarly, retrievals from McBride et al. (2011) were based on hyperspectral 
zenith radiances, and compared to two MODIS overpasses for overcast water clouds at the 
ARM Oklahoma site in 2007. 
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Fig. 2. AERONET cloud-mode site locations (red squares). Among these cloud-mode sites, 
31 providing zenith radiance data at 1640 nm wavelength are listed in Table 1 . 
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Cloud optical depth 


Fig. 3. Dependencies of normalized zenith radiance on cloud optical depth at 870 and 1640 nm 
wavelengths for cloud effective droplet size of 8 and 16 pm. Surface albedo values are 0.3 at 
870 nm, and 0.25 at 1 640 nm. Solar zenith angle is 45°. 
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Cloud opticol depth 



X-km X-km 


Fig. 4. (a) A Large Eddy Simulation-generated cloud optical depth field used to calculate zenith 
radiances that are then in turn used to test the retrieval method. Taking the transect along the 
thick black-dashed line (around 3 km in the y direction) gives the distribution of effective radii in 
the x- and vertical directions in (b). 
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Fig. 5. Intercomparison between the true and the retrieved liquid water paths using box plot (a) 
and scatter plot (b), for the simulation in Fig. 4. The bottom and top of each box represent the 
25% and 75% quartiles, and the line inside the box represents the median (Tukey, 1977, p. 
41-43). The whiskers mark the “accepted range,” which is defined as 1 .5 times the interquartile 
distance. The open circles outside the accepted range are considered outliers. The error bars 
in (b) represent one standard deviation of retrievals from the 40 realizations in the retrieval 
process, while the co-plotted dashed lines represent a typical uncertainty of 30 gm -2 in liquid 
water path retrieved from microwave radiometer measurements, (c) and (d) are the same as 
(a) and (b), respectively, but for cloud optical depth. 
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Fig. 6. Scatter plots of retrieved versus true effective radius from simulations. Two “true” effec- 
tive radii are calculated: T reff constLWC from the assumption of constant liquid water content in the 
vertical; 7 eff from averaging droplet radii over all levels in the cloud layer. In addition, two sets 
of retrieved effective radii are used: one is based on simulation at 201 m resolution, the other 
is at 67 m resolution, (a) and (b) are retrievals at 201 m resolution plotted against T reff constLWC 
and 7 eff . The error bars represent the standard errors (discussed in Sect. 2). (c) Box plots for 
the two types of the true effective radii, and the two sets of retrievals. 
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Fig. 7. Box plots of the effective radii at locations along the transect at 3.1 km in the Y direc- 
tion, as shown in Fig. 4. (a-d) are one continuous plot broken into four panels. Dots are the 
corresponding retrievals. 
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Fig. 8. (a) Scatter plot of liquid water paths retrieved from the ARM microwave radiometer 
(MWR; Turner et al. , 2007) versus coincident cloud-mode measurements at 1.5 min temporal 
resolution for all liquid water clouds during May 2007-June 2008. Cloud-mode LWP values are 
calculated from retrieved cloud optical depth and effective radius, assuming that cloud droplet 
concentration is approximately constant and that liquid water content increases linearly with 
height, (b) The corresponding occurrence counts using a bin size of 25gm“ 2 . 
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Cloud effective radius (nm) 

Fig. 9. Intercomparison of cloud effective radius retrieved from cloud mode to those from the 
ARM Mace dataset by showing (a) the histograms of the occurrence counts, and (b) the box 
plots. Note that cloud-mode retrievals represent a 1 .5 min average droplet radius, while the 
ARM Mace retrievals represent a 5 min average radius. 
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Fig. 10. (a) Attenuated backscatter signal from micropulse lidar on 15th June 2007. A time 
series of effective radius (b) and cloud optical depth (c) from cloud-mode measurements and 
the ARM Min dataset. The error bars represent one standard deviation. In (b), we omit the 
default values of 8 pm reported in the ARM Min retrievals. In (c), cloud optical depth from the 
ARM Min dataset was truncated at 100 for plotting purposes. 
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Fig. 11. Scatter plots of 1 .5-min average cloud optical depth (a) and effective radii (b) retrieved 
from cloud mode versus those from the ARM Min dataset. Error bars for the ARM Min dataset 
represent one standard deviation of retrievals during 1 h time windows. 
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Overpass index 


Fig. 12. Box plots of cloud effective droplet radii retrieved from Terra/Aqua MODIS in 2007. 
Coplotted dots represent cloud-mode 1.5 min average effective radius during the 1 h time win- 
dow centred at the MODIS overpass time. Detailed information for each overpass can be found 
in Table 5. 
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Fig. 13. Same as Fig. 1 , but for intercomparison to MODIS-retrieved cloud effective droplet radii 
only. Cloud-mode retrievals are based on eight cases listed in Table 5. 
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Fig. A1. Plots of the retrievals versus the truth, assuming that the retrievals agree to the truth 
within an error of (a) 0.5 gm, (b) 1.0 gm, (c) 1.5gm, (d) 2.0 gm. The corresponding correlation 
coefficient (Ft) is listed in each plot. 
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